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Abstract 


We report the detailed history of spin-period changes in five intermediate polars (DQ Herculis, AO Piscium, FO 
Aquarii, V1223 Sagittarii, and BG Canis Minoris) during the 30—60 yr since their original discovery. Most are 
slowly spinning up, although there are sometimes years-long episodes of spin-down. This is supportive of the idea 
that the underlying magnetic white dwarfs are near spin equilibrium. In addition to the ~40 stars sharing many 
properties and defined by their strong, pulsed X-ray emission, there are a few rotating much faster (P < 80 s), 
whose membership in the class is still in doubt—and who are overdue for closer study. 


Unified Astronomy Thesaurus concepts: Cataclysmic variable stars (203); Classical novae (251); Close binary stars 
(254); DQ Herculis stars (407); Interacting binary stars (801); Novae (1127); Stellar accretion (1578); Stellar 


accretion disks (1579) 


1. Introduction 


Intermediate polars (IPs, also called DQ Her stars after the 
prototype) are magnetic cataclysmic variables with stable 
periodic signals in optical and X-ray light, and periods typically 


31 Visiting Astronomer, Cerro Tololo Inter-American Observatory, National 
Optical Astronomy Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., (AURA) under cooperative 
agreement with the National Science Foundation. 
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in the range of 1-30 minutes. These periods come from rotation 
of the radially accreting, magnetic white dwarf (WD), although 
some also show “sideband” signals which arise from interaction 
between the spin and orbital clocks. The first was found long ago, 
in the remnant of Nova Herculis 1934 (Walker 1956, 1961). 
X-ray telescopes since 1980 have revealed many more, because 
most IPs radiate most of their energy in X-rays, as accreting 
matter plunges radially to the WD surface. Patterson (1994, 
hereafter P94) reviews these stars, and Table 1 of Norton et al. 
(2004, hereafter NWS) presents a more complete list of class 
members. The NASA website created by Koji Mukai?” contains 


32 The Intermediate Polars, https: //asd.gsfc.nasa.gov /Koji.Mukai /iphome/ 
iphome.html. 
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by far the most up-to-date and useful online list of the ~50 
class members and their individual properties. Much of the 
material in Mukai’s (2017) review of X-ray emission in 
cataclysmic variables is also very pertinent to IPs. 

By tracking period changes from year to year, one can in 
principle measure torques on the rotating WD. This can 
constrain the accretion rate and the WD magnetic moment 
(P94; NWS). Several authors have attempted this, based on 
time-series photometry obtained over a baseline of several 
years. This has been sufficient to yield a rough estimate for the 
several stars studied: they change their pulse periods on 
timescales [=P/(dP/dt)] around 10° yr. 

We have carried out such programs over many years, most 
recently with the globally distributed small telescopes of the 
Center for Backyard Astrophysics (CBA; Patterson et al. 2013). 
Our baselines are very long, usually from the discovery year to 
the present. About 30 of the ~50 known class members are in 
our archives: nearly all IPs brighter than V ~ 17. Most are 
monitored 3-15 times per year, with special care to obtain 
timings early and late in each observing season (which 
eliminates errors in counting cycles between consecutive 
seasons). 

Faithful tracking of known IPs appears to be less glamorous 
than the discovery of new class members. Many of the 
published studies base their period estimates on just one 
observing season—sufficient to establish the stability of the fast 
signal, but not to measure period change. In the course of this 
work, we have also found that some of the published spin 
ephemerides spanning more than ~3 yr are incorrect. The main 
reason is cycle-count errors between years, because observers 
tend to disfavor the poorer observing conditions of the early 
and late season. In addition to frustrating the search for dP/dt, 
lack of a reliable long-term ephemeris, or at least accurate 
period estimation, also hampers interpretation of data at other 
wavelengths. But full publication of our results will take years 
to complete, because of volume (>9000 nights so far), and 
rapid discovery of new IPs makes our task ever more daunting. 

So we here present a summary of the period history of five IPs 
with the longest baseline of observation. These typify the patterns 
and timescales found in other class members, but are more clearly 
defined, because the baseline is longer (at least 35 yr). 


2. Individual Stars 
2.1. Measurements 


In studies of period change, it is standard practice to present 
“O-C curves,” which represent a star’s departure from a 
strict constant-period ephemeris (e.g., Breus et al. 2019, or 
Kreiner 1971 for a very extensive application to close binaries 
generally). For orbital period changes, it is essentially 
mandatory, because the period change is so small as to require 
the full decades-long baseline to yield a tiny measurable effect. 
Such studies always present precise timings of individual 
events (e.g., maximum/minimum light, mid-eclipse) and then 
fit timings with a constant period or low-order polynomial (P, 
P, maybe Ë). 

This technique flounders when the event being timed is not 
precisely defined (“maximum light’), or is not observed 
sufficiently often to establish cycle count with certainty, or 
when the period changes too fast. Intermediate-polar spin 
history brings each of these problems into play. In addition, as 
this study will show, the spin periods of these IPs typically 
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wander on timescales of years, for no clear reason. Thus the È 
and Ë terms of a polynomial fit may be mere accidents of the 
observing interval. Finally, because our study involves 
contributions from ~25 different telescopes over ~40 yr, the 
merging of data on a common scale can be a problem. 

For these reasons, we present the data not as O-C curves, but 
as period versus time [P(t)]. For each star in our program, the 
original data consists of nightly time-series photometry with 
typically 20s integrations and 3-8hr in length. We first 
calculate the average nightly power spectrum of the best and 
longest “nights” (usually 8—20hr long). The frequencies 
revealed will typically be the orbit Q, the spin pulse w, its 
lower orbital sideband w—Q, and several higher harmonics. This 
is always low-noise (because of the averaging) but low 
resolution. We then look for the most densely sampled long 
segments (usually 10-50 nights long) and calculate the power 
spectrum for each. These appear to resolve, with high accuracy, 
all the fine-structure in the power spectrum: the main signals w 
and Q, plus many weak signals which take the form nw—mQ, 
where m = 0, 1, ..., n (similar but not identical to the fine- 
structure of “superhump” signals, as discussed by Skillman 
et al. 1999).*° The spread in terrestrial longitude insures against 
selecting the wrong daily alias frequency. In many years we 
take care to obtain observations early and late in the ~6 month 
observing season; this insures against selecting the wrong 
yearly alias. 

We then return to the individual nights, or consecutive-night 
segments, and do a synchronous summation at the dominant 
high-frequency signal (usually w, but sometimes w—Q) to yield 
times of maximum light. Then we collect all the maximum- 
light timings (typically 15-50) over 3 yr intervals (average 
baseline ~2.3 yr), make a weighted linear fit, and record the 
periods as running averages over those baselines. This appears 
to give good accuracy. For example, a 20 minute signal will 
execute ~60,000 cycles during 3 yr (i.e., an interval of 
~2.3 yr); if each timing is accurate to 0.07 cycles (a reasonable 
but conservative estimate), then the period is measured to an 
accuracy of 0.0013 s. With many (>15) such timings, the error 
shrinks to ~0.0008 s. These numbers are typical for most of 
our data. The running-average periods are shown with their 
yearly midpoints in Table 1. 

To illustrate this point, Figure 1 presents a traditional O-C 
diagram for the 913 s signal in BG Canis Minoris, which has 
previously been flagged as an IP with erratic behavior (Kim 
et al. 2005; Bonnardeau 2016). The observations span 38 yr 
with no uncertainty in cycle count. Generally, the “sheds 
water” shape of the curve shows that the period decreases 
throughout. But the rate of decrease is not constant, and there is 
no simple mathematical expression to describe it. P(t), as 
described, is probably a better way to render a long history; 
compare the O-C curve in Figure | with the corresponding P(t) 
representation’ for BG CMi in Figure 2. 

Finally, our identification of w as the true spin frequency is 
not always guaranteed to be correct. Only X-ray time series can 
unambiguously distinguish w from w—Q; and even then, there 





33 The commonality of the -Q feature(s) presumably arises not from any deep 
similarity, but just from the fact that both IPs and common superhumps rotate 
prograde with respect to the orbit. 

4 The pulse-period history of the high-luminosity pulsating X-ray sources 
(“X-ray pulsars”) is always represented this way (e.g., the many figures and 
superb analysis in Bildsten et al. 1997). The reason is that the periods change 
much too fast for O-C analysis; neutron stars are easy to spin up (or down)! 
With infrequent observation, the same applies to IPs. 
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Table 1 
History of Spin Periods 
DQ Her AO Psc FO Agr V1223 Sgr BG CMi 
Period (s) Period (s) Period (s) Period (s) Period (s) 
Year (71+) Year (858+) Year (1254+) Year (794+) Year (913+) 
1955.48 0.065858 1980.5 0.6893 1982.5 0.4487 1981.4 0.3803 1982.7 0.5055 
1957.06 0.065804 1981.5 0.6860 1983.9 0.4474 1982.8 0.3808 1983.6 0.5041 
1958.20 0.065796 1983.5 0.6838 1985.5 0.4530 1984.1 0.3827 1984.4 0.5017 
1959.1 0.065774 1986.0 0.675 1986.7 0.4537 1985.7 0.3832 1985.6 0.4978 
1968.07 0.065564 1988.1 0.671 1987.7 0.4514 1988.6 0.3863 1987.1 0.496 
1969.65 0.065556 1990.2 0.6670 1989.0 0.4522 1998.5 0.3934 1988.6 0.492 
1970.2 0.065540 1993.3 0.661 1990.4 0.4505 2001.0 0.3959 1989.6 0.4895 
1971.68 0.06550 1996.3 0.6550 1991.4 0.4494 2004.0 0.3983 1990.2 0.4877 
1973.25 0.065456 1998.3 0.6530 1992.2 0.4480 2005.7 0.3986 1991.4 0.4867 
1975.06 0.065435 2000.3 0.6517 1993.4 0.4459 2007.2 0.4023 1993.7 0.4802 
1976.27 0.065413 2001.3 0.6500 1995.4 0.4394 2008.5 0.4018 1995.4 0.4756 
1977.1 0.06540 2002.3 0.6482 1996.2 0.4334 2010.0 0.4029 1996.9 0.4723 
1978.30 0.065374 2003.3 0.6474 1997.2 0.4272 2011.3 0.4043 1998.8 0.4693 
1980.77 0.065342 2004.3 0.6446 1998.2 0.4206 2012.5 0.4045 2000.4 0.4707 
1982.80 0.065313 2005.3 0.6428 1999.4 0.4138 2013.7 0.4052 2001.5 0.4701 
1985.5 0.065268 2006.3 0.6404 2000.2 0.4066 2014.5 0.4054 2002.6 0.4717 
1989.5 0.065203 2007.2 0.6384 2001.2 0.3982 2016.2 0.4066 2003.6 0.4697 
1992.0 0.065170 2008.3 0.6380 2002.2 0.3916 2017.2 0.4069 2005.3 0.4717 
1993.5 0.065158 2009.3 0.6358 2003.2 0.3841 vee tee 2006.7 0.4712 
2000.70 0.065050 2010.3 0.6350 2004.2 0.3820 tee tee 2008.5 0.4702 
2001.3 0.065037 2011.3 0.6321 2005.2 0.3753 te vee 2009.7 0.4703 
2002.3 0.065028 2012.3 0.6308 2006.2 0.3714 vee vee 2010.5 0.4681 
2007.44 0.064972 2013.3 0.6285 2007.2 0.3614 nee ves 2011.4 0.467 
2008.34 0.064953 2014.3 0.6285 2008.2 0.3560 ve vee 2013.7 0.4658 
2009.58 0.064925 2015.3 0.6256 2012.2 0.3396 tee tee 2016.3 0.4653 
2010.59 0.064915 2016.3 0.6239 2014.2 0.3311 nee ves 2017.4 0.4644 
2011.94 0.064896 2017.3 0.6219 2015.2 0.3324 
2012.85 0.064888 tee ves 2016.2 0.3360 
2013.52 0.064883 vee tee 2017.2 0.3379 
2014.65 0.064864 vee tee 
2017.0 0.064832 
5.0 To > 2.2. DQ Herculis (Nova Herculis 1934) 
[ BGCMI ] DQ Her was regarded as “the nova of the century” for most 
L ate 41982-2018 4 of the twentieth century. Among its several first-ever contribu- 
a 2.5 :° Pe 3 a tions to cataclysmic-variable science was the discovery of 
a - Pi “ee 4 strictly coherent 71 s pulses in the light curve (Walker 1956, 
D | r E "ey | 1961). This has been tracked continuously ever since: the most 
O OF %» = recent study is that of Wood et al. (2005). The P(t) history is 
(®) [ g my | given in Table 1, and shown in Figure 2. The period has 
ò T 2 *. ] decreased continuously since 1954, with a rate of period 
25+ ¢ s decrease declining from 26 to 12 us yr’. 
E Lie S] DQ Her is exempt from one of the caveats expressed in 
- a | Section 2.1, because the cycle count is known with certainty 
50 bee ag = ee pe a a from 1954 to the present. On the other hand, the timing of the 
0 0.5x10° 4.0x10° 71 s signal is known to present a large, systematic dependence 
on orbital phase (Warner et al. 1972; Patterson et al. 1978). So 
E (CYCLES) 


Figure 1. O-C diagram for maximum light of the 913 s signal in BG CMi, 
relative to an assumed period of 913.48 s. The general shape indicates a spin-up 
(period decrease) over the 38 yr of observation, but with a rate that is not 
consistent with a simple polynomial fit. 


can be a 2x ambiguity for a two-pole accretor. In our time- 
series photometry, we track the wanderings of the usually 
dominant high-frequency optical signal. Other work is 
generally needed to reveal if it is actually w, or w-Q, or 2w, 
or (for an exotic geometry) some other combination. 


a reliable timing should exclude the eclipse and be averaged 
over the remainder of the orbit. This is true for most of our 
data; and since each period estimate is based on at least a dozen 
(usually >20) nights of observation, this concern should not 
affect our result. 

The accuracy of the DQ Her timings is much greater than 
that of the other stars considered here. This is because the 
signal is ~13x faster, because there are no sideband signals to 
confuse matters, because the signal is accurately sinusoidal, 
and because the observing season is long. We estimate that 
each 3 yr period estimate is accurate to ~0.003 ms. 
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Figure 2. Period vs. time for the five IPs with longest duration of observation. The baselines for most points are in the range of 500-800 days, and the resultant errors 
(which represent an estimated phase uncertainty of ~0.07 cycles over those baselines) are about the size of the symbols. 


2.3. AO Piscium (H2252-035) 


Following DQ Her, AO Psc was the second linchpin in the 
discovery of IPs. Griffiths et al. (1980) found the star to be a 
strong X-ray-emitter in the HEAO A-3 data. Subsequent study 
showed stable optical signals at 3.6 hr and 859s (Patterson & 
Price 1981), and an X-ray signal at 805s (White & 
Marshall 1981). These studies identified 805s as the WD’s 
spin period, 3.6 hr as the orbital period, and 859 s as the lower 
(in frequency) orbital sideband of the spin frequency. Both of 
the fast periods are actually present in the optical photometry, 
and their beat frequency equals the orbital frequency to within 
1 part in 10°. This three-period structure came to be the 


standard for IPs, although one of the fast periods is sometimes 
missing (within observational limits). 

In recent years it has been possible to measure the X-ray 
spectrum to good precision, and IPs are generally found to be 
very hard sources, with temperatures exceeding 40 keV. This 
agrees well with the theory that the main power source is radial 
accretion from gas accreting along magnetic-field lines, high 
above the magnetic pole. The X-ray frequency fy is taken to be 
the true spin frequency. For the well-studied cases, the 
dominant fast optical signal occurs either at fy or at the lower 
orbital sideband fy — for», as expected for prograde rotation of 
the WD (i.e., no certifiably retrograde rotators are known). 
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The spin-period history of AO Psc’s dominant optical (859 s) 
signal is tabulated and tracked in Table | and Figure 2. For the 
last 40 yr, the accreting WD has been spinning up at a rate of 


~1.7 msyr’. 


2.4. FO Aquarii (H2215—086), V1223 Sagittarii (4U185 1-31), 
and BG Canis Minoris (3A0729+103) 


After the discovery of AO Psc, three other hard X-ray 
sources were quickly found to coincide with stars showing 
spectra and light curves characteristic of cataclysmic variables, 
and with similar fast periods found in optical photometry. 
Steiner et al. (1981) found a 794s period in V1223 Sgr, 
Patterson & Steiner (1983, see also Shafter & Targan 1982) 
found a 1254s period in FO Agqr, and McHardy et al. (1984) 
found a 913s period in BG CMi. Several papers have tracked 
period changes since then (V1223 Sgr: Jablonski & Steiner 
1987; FO Aqr: Patterson et al. 1998, Littlefield et al. 2019; BG 
CMi: Patterson & Thomas 1993; see also Norton et al. 1992). 
After collecting the published data and adding our own ~400 
timings, we find the P(t) behavior recorded in Table | and 
shown in Figure 2. 

Figure 2 is a representative collection. Most IPs decrease 
their periods, as one might expect, since they are probably 
accreting gas from a disk, which has higher specific angular 
momentum than the WD. But some show episodes of period 
increase, which can last as long as 35 yr or more (V1223 Sgr). 
Since the observed timescales of period change (P/P from 
Figure 2) are typically near 10° yr—far shorter than the lifetime 
of the stars—it is generally thought that the stars are near “spin 
equilibrium,” where they vacillate between episodes of spin-up 
and spin-down. That is very likely true, although it has never 
been proven. If they do vacillate, it can be on timescales as long 
as 50 yr. 


3. Spin-period Change in Theory 


The prevailing theory for spin-period change in accreting, 
magnetic compact stars is that of Ghosh & Lamb (1979). This 
is reviewed and applied to IPs by NWS, P94, Lamb & 
Patterson (1983), and Mukai (2017). Stars with decreasing 
periods are deemed to be “slow rotators,” with the spin-up 
matter torque of accreting gas exceeding the spin-down torque 
of the WD’s magnetic-field lines entangled in the outer, slowly 
rotating disk. But many consecutive years of spin-up will move 
the star to “fast rotator” status, where the WD’s field lines 
entangle in the outer disk and slow the star down. Thus is 
created a spin equilibrium, although modulated by any changes 
in mass-transfer rate—endemic to all cataclysmic variables. 

It should be stressed that the dP/dt values (the slopes) in 
Figure 2 do not, in this theory and any plausible theory, 
represent actual evolution times—but merely some accidental 
feature of the current era (possibly mass-transfer rate; the WD 
magnetic moment is also critical, but is assumed constant). 

Fast-rotator status tends to inhibit accretion, since it invokes 
a centrifugal barrier. Therefore we expect that stars will be 
fainter during the fast-rotator phase—and thus predict that most 
known IPs should be in their slow-rotator phase (spinning up). 
As apparently observed. Figure 16 of P94 shows the 
general idea. 
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4. V1223 Sagittarii and Spin-down Episodes 


In Figure 2, V1223 Sgr seems to be an exception. All the 
other stars show period decrease, with at most small and short- 
lived episodes of period increase. And this appears to be 
generally true for the other several dozen stars in our program 
(with more fragmentary data; there may be some interesting 
exceptions). Now our general idea is that these stars, 
supposedly near spin equilibrium, spend comparable times 
spinning up and down; so why is this a surprise? Because a 
WD spinning down should be in a state of low luminosity, as 
its flailing magnetic-field lines drag in the slowly rotating disk. 
This explains why most IPs are spinning up (because the stars 
are harder to detect in a low-luminosity state). Yet V1223 Sgr 
appears to be one of the more luminous IPs (Beuermann et al. 
2004; Schwope 2018). 

So there are still some mysteries to unpack on this subject. 
One possibility lies in the star’s long-term history of brightness 
fluctuations, which seem to be much larger than those certified 
in other IPs. Garnavich & Szkody (1988) report a long-lasting 
“low state” in the 1940s, and Simon (2015, Figures 1 and 3) 
shows that this large variability in brightness is quite 
characteristic of the star. 

Littlefield et al. (2016, 2019) have studied the history of such 
period changes in FO Aqr, and found evidence for the expected 
correlation between spin-down and low states. This certainly 
has promise, although will be hampered by erratic variability, 
any looseness in the classification of “low states,” and the long 
timescale required for period changes to be manifest. 


5. Short-period Cousins? 


Among the ~50 known IPs, ~35 have been X-ray-selected. 
They all have common properties: strong and hard X-rays, 
high-excitation emission lines, spin periods exceeding 3 
minutes, and the presence of sideband signals. The common- 
alities are sufficiently extensive to warrant lumping into one 
class: IPs. 

But there are a few other cataclysmic variables with stable, 
very short periods: WZ Sge at 27.87s; AE Aqr at 33.08 s, 
V533 Her at 63.63 s, V455 And at 67.62 s; and possibly DQ 
Her itself at 71.06 s. They each have unique quirks which are 
the subject of many research papers, and it remains unknown’? 
if their underlying physics is predominantly that of the IPs. 
Their principal disqualifier is the absence of strong, hard, 
pulsed X-rays. But on the other hand, that may merely be the 
result of fast rotation (see Section 5.7 of Mukai 2017). 
Considering that most of these fast periodic signals were 
discovered in the 1970s—before any of the well-credentialed 
IPs were found—it seems likely that a well-designed search for 
more candidates would be fruitful. 


6. Are Nova Eruptions Relevant? 


It should be remembered that IPs have an extra piece of 
physics (thermonuclear eruptions) not available to their 
accreting neutron-star cousins. This could entail some other 
sources of spin-up: either an accretion torque specifically 


33 DQ Her itself is generally included in the IP class, on grounds of seniority, 
period stability, pulsed high-excitation emission lines, and a good excuse for 
concealing X-rays (a binary inclination close to 90°). The others are usually not 
included in lists—though possibly because lists are often prepared by X-ray 
astronomers. Without that energy bias, it is possible that all but WZ Sge would 
qualify. 
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associated with the elevated M following a nova outburst 
(which would make IPs temporarily “slow rotators” and thus 
spin up as their magnetospheres are squashed), or a slow 
contraction (conserving angular momentum) of a WD heated 
by a recent outburst. Another possibility is that some of the 
large outflow from the WD during the main outburst, or 
immediate aftermath, proceeds along magnetic-field lines and 
thereby slows the WD rotation. These could be called “Blame 
It On The Bossa Nova” theories (Gormé et al. 1963). The 
energies involved seem reasonable: a century of spin-up at 
1 ms yr | taps only 0.01% of the total energy (10*° erg) of the 
nova explosion. And we note from Figure 2 that the observed 
P(t) in DQ Her, interpreted as exponential decay, appears to 
suggest a time constant of ~60 yr—roughly the age of the 
postnova. 


7. Summary and a Look Ahead 


1. Period-versus-time (Figure 2 and Table 1) is for most 
purposes the best way to illustrate the period changes— 
rather than the more traditional O-C diagram, which is 
hard to interpret when the changes are not monotonic. It 
should be useful to observers attempting to phase their 
data, and to theorists trying to understand accretion 
torques in these and related stars (e.g., X-ray pulsars). We 
would be happy to reduce our analysis burden by 
furnishing data on any of the ~30 stars to interested 
researchers. 

2. For four of five stars reported here, and for most IPs with 
shorter baselines of observation, spin-up is the general 
rule, and the shortness of the timescale (P/P ~ 10° yr) 
suggests that the WDs are near spin equilibrium. But that 
requires episodes of spin-down, and the causes for such 
episodes are not yet known. An extensive and calibrated 
record of visual or X-ray brightness, over a baseline of 
years and supplemented by period data of the type 
reported here, may reveal those causes. However, the spin 
history of accreting neutron stars (Bildsten et al. 1997) 
suggests that the relationship between dP/dt and bright- 
ness may prove to be quite complex (e.g., Perna et al. 
2006). 

3. In the present century, almost nothing new has been 
learned about the fastest IPs. This signifies not the 
maturity of the subject, but mainly the lack of human 
attention. The stars themselves seem cooperative: some 
are very nearby, some shout for attention via classical- 
nova outbursts, and some have a long history of previous 
work which has never been well-digested in the context 
of what is now known about the slower IPs. That might 
be a fruitful subject area for the 2020s. 


To span 40 yr, you have to keep observing for 40 yr. In 
addition to all the small telescopes contributing to this paper, 
there were also ~500 nights on meter-class telescopes, mainly 
from Kitt Peak, Cerro Tololo, and Lick Observatory. We owe a 
great debt to those mountain staffs for all the technical (and 
medical!) help they provided. Jules Halpern, Ivan Andronov, 
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and especially the anonymous referee have helped us improve 
the paper. On the financial side, grants from the NSF (AST- 
1615456 and AST-1908582), the Research Corporation, and 
the Mount Cuba Astronomical Foundation have been key 
elements in this work. The work also includes observations 
obtained with the Mittelman Observatories 0.5 m telescope at 
New Mexico Skies in Mayhill, New Mexico. Additionally, we 
would like to thank the Mittelman Family Foundation for its 
generous support of and substantial impact on astronomy at 
Middlebury College. 
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